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(vq ■ ABSTRACT 

Aims. Abundances of Mg isotopes 24 Mg, 25 Mg, and 26 Mg can be used to test models of chemical enrichment of 
St ■ interstellar/intergalactic gas clouds. Additionally, since the position of the Mgn AA2796, 2803 A lines is often taken as 

a reference in computations of possible changes of the fine-structure constant a, it should be clarified to what extent 
l^-} ' these lines are affected by isotopic shifts. 

Methods. We use a high-resolution spectrum (pixel size ~1.3 km s _1 ) of the quasar HE000 1-2340 observed with 
the UVES/VLT to measure Mg isotope abundances in the intervening absorption- line systems at high redshifts. Line 
profiles are prepared accounting for possible shifts between the individual exposures. In the line fitting procedure, the 
lines of each ion are treated independently. Due to unique composition of the selected systems - the presence of several 
transitions of the same ion - we can test the local accuracy of the wavelength scale calibration which is the main source 
of errors in the sub-pixel line position measurements. 
Results. In the system at z a b s = 0.45 which is probably a fragment of the outflow caused by SN la explosion of 
high-metallicity white dwarf(s) we measured velocity shifts of Mgn and Mgl lines relative to other lines (Fel, Fell, 
Cal, Can): AV Mg ii = —0.44 ± 0.05 km s _1 , and AVmsi = —0.17 ± 0.17 km s -1 . This translates into the isotopic 
ratio 24 Mg: 25 Mg: 26 Mg = (19 ± 11) : (22 ± 13) : (59 ± 6) with a strong relative overabundance of heavy Mg isotopes, 
( 25 Mg+ 26 Mg)/ 24 Mg = 4, as compared to the solar ratio 24 Mg: 25 Mg: 26 Mg = 79:10:11, and ( 25 Mg+ 26 Mg)/ 24 Mg = 0.3. 
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In the systems at z a b s = 1-58 and z a b s = 1-65 enriched by AGB-stars we find only upper limits on the content of heavy 
Mg isotopes ( 25 Mg+ 26 Mg)/ 24 Mg < 0.7 and ( 25 Mg+ 2 %g)/ 24 Mg < 2.6, respectively. At z abs = 1.58, we also put a 
strong constraint on a putative variation of a: Aa/a — ( — 1.5 ± 2.6) x 10~ 6 which is one of the most stringent limits 
obtained from optical spectra of QSOs. We reveal that the wavelength calibration in the range above 7500 A is subject 
to systematic wavelength-dependent drifts. 



Key words, line: profiles - methods: data analysis - galaxies: abundances - quasars: absorption lines - quasars: indi- 
vidual: HE0001-2340 



£^ ■ 1. Introduction sions of CO white dwarfs where 22 Ne is the third most 

^— ' ' abundant nuclide (Jose et al. 1999, 2004), or in outflows 

Magnesium is one of a few elements for which the isotope from AGB-stars where the so-called hot-bottom burning 



abundances can be measured from astronomical spectra. operates (Karakas et al. 2003, 2006, 2010). Direct measure- 
It has three stable isotopes - the alpha-nucleus Mg and ments of Mg isotopic ratio in different objects allow us to 
the neutron-rich nuclei b Mg and b Mg. In the solar pho- ide ntify the sources of the chemical enrichment and thus to 
tosphere, the magnesium isotopes are mixed m the proper- underst and the chemical evolution of these objects. 



X 

^ ' tion 24 Mg: 25 Mg: 26 Mg = 78.99:10.00:11.01 (Morton 2003) 



High-resolution spectral observations of other stars in the Beyond the study of the chemical enrichment of the 

Milky Way show that the content of heavy isotopes can intergalactic gas, measurements of the Mg isotopic abun- 

be higher, e.g., 24 Mg: 25 Mg: 26 Mg = 48:13:39 (Yong et al. dances in the intervening clouds at high redshifts are closely 

2003, 2006). An enhanced content of 25 Mg, 26 Mg - up to related to the current debate on a putative variation of 

24 Mg: 25 Mg: 26 Mg = 60:20:20 - was also measured in some the fine-structure constant a (for a review, see, e.g., Uzan 

presolar spinel grains (Zinner et al. 2005; Gyngard et al. 2010). Strong resonance lines of the Mgn doublet (AA2796, 

2010) which are thought to be related to the interstellar 2803 A) are often used as a reference ('anchor') relative to 

dust. 25 Mg and 26 Mg are produced from 24 Mg via proton which the positions of all other ions and, hence, the value 

capture in the MgAl chain and from 22 Ne via alpha-capture of Aa/a = (a spaco - aq a b)/ai a b is calculated (Murphy et 

and are supposed to be enhanced in ejecta of nova explo- al. 2001; Chand et al. 2004). However, if Mg isotope abun- 

dances in quasar absorbers differ from the solar value, this 

* Based on observations performed at the VLT Kueyen tele- manifests itself in a shift of the reference lines and, hence, 

scope (ESO, Paranal, Chile), the ESO programme No. 183. A- can mimic variations of a (Levshakov 1994; Ashenfelter et 

0733 al. 2004). 
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In our recent paper (Levshakov et al. 2009, hereafter 
Paper I) we reported on a tentative detection of a blueward 
shift of the Mgn AA2796, 2803 A lines in several metal- 
rich absorbers at z a b s ~ 1.8, which may be caused by an 
enhanced content of 25 Mg and 26 Mg in the absorbing gas. 
However, this was a rather qualitative result since the avail- 
able quasar spectra had neither a sufficiently high spectral 
resolution nor a calibration accuracy needed to evaluate a 
particular isotopic abundance ratio. 

Heavy Mg isotopes have shorter wavelengths than 24 Mg, 
and an enhanced content of 25 Mg and 26 Mg in the mixture 
24 Mg+ 25 Mg-|- 26 Mg can be revealed from a blueward shift 
of the barycenter of Mg lines relative to its value based 
on the laboratory measurements (i.e., with the solar abun- 
dance ratio). The maximum possible isotope shift for the 
Mgn lines is —0.8 km s _1 - when all absorbing Mg con- 
sists of 26 Mg. For the isotope ratios measured in the Milky 
Way giants and in presolar spinel grains the offset of Mg II 
lines would be several times smaller - from —0.25 to —0.1 
km s _1 . Since the limiting pixel size which can be achieved 
in spectral observations of such faint sources as quasars is 
of ~ 1 km s _1 , all expected blueward shifts of the Mgn 
lines caused by the presence of neutron-rich species are at 
the sub-pixel scale. To detect such shifts we need a very 
bright object to reach a high S/N ratio ( > 30) even at 
high spectral resolutions. Besides, in order to distinguish 
isotope shifts of the Mg n lines from kinematic shifts caused 
by the velocity-density inhomogeneities inside the absorber, 
a special attention should be paid for the selection of ab- 
sorbers: appropriate systems ought to contain narrow lines 
of several low-ionization metal ions with simple symmetric 
profiles. 

While performing the measurements on the sub- 
pixel level, the calibration accuracy becomes of the ut- 
most importance. Recent tests with iodine cells for both 
Keck/HIRES and VLT/UVES spectra show that calibra- 
tion uncertainties can result in wavelength offsets up to 
1000 m s" 1 (Griest et al. 2010; Whitmore et al. 2010). 
The local wavelength calibration can be verified by com- 
parison of different lines of the same ion like, e.g., Sin 
AA1260, 1304, 1526 A, or Fen AA2344, 2382, 2600 A. This 
again requires the presence of such absorption systems in 
quasar spectra. 

Accounting for all these conditions, a bright quasar 
HE0001-2340 (V = 16.6, z cm = 2.28) discovered in course 
of the Hamburg/ESO survey (Wisotzki et al. 1996; Reimers 
et al. 1996; Wisotzki et al. 2000) seems to be a perfect 
target: several systems with strong and narrow Mgn lines 
accompanied by lines of other low-ionization ions arc identi- 
fied in its spectrum. These are the systems at z a b s = 0.45207 
(D'Odorico 2007), z abs = 1.5864 (Chand et al. 2004), z abs 
= 1.6515 (Paper I) which are suitable for measuring the 
Mg isotope abundances, and a sub-DLA system at z a b s = 
2.1871 (Richter et al. 2005) exhibiting several sets of ab- 
sorption lines of the same ion which can be used to test the 
local wavelength calibration. 

In this paper the subsequent Sect. [2] presents observa- 
tions and data reduction, the analysis of the line profiles is 
described in Sect. [31 the results obtained are discussed in 
Sect. HJ and Sect. [5] sums up the present study. 



2. Observations 

The observations of the quasar HE0001-2340 were acquired 
with the UV- Visual Echelle Spectrograph (UVES) at the 
VLT 8.2-m telescope at Paranal (Chile) on 5 nights in 
September 20-24, 20093. The journal of observations is re- 
ported in Table [T] together with additional information. 
Some exposures were taken with the standard dichroic 
beam splitter DIC2 (settings 437+760 and 420+700), and 
some with the DIC1 (settings 390+580) thus providing the 
total wavelength coverage between 350 and 860 nm. The 
wavelength ranges covered by individual settings are given 
in Tabled The slit width was set to 0.7 arcsec for all obser- 
vations providing a resolving power of f» 65554 + 3868. This 
slit width is narrower than usually used for observing high- 
redshift quasars. During the observations the seeing was be- 
tween 0.5 arcsec to 1 arcsec as measured by the Differential 
Image Motion Monitor (Sarazin & Roddier 1990) but gen- 
erally was slightly better at the telescope. The major differ- 
ence with all previous observations was that the CCD pixels 
were read without binning. Data reduction was performed 
with the last version 4.4.8 of the UVES pipeline (Larson 
& Modigliani 2009). The optimal extraction method was 
adopted to extract the flux in the cross dispersion direc- 
tion. The 2D images of the long-slit calibration lamp and 
long-slit flat-field have been bias subtracted but not flat 
fielded (in the UVES, a pinhole lamp is used for the loca- 
tion of the orders which are curved and somewhat tilted 
upward). In the final spectrum the pixel size is of « 1.25 
km s- 1 (FWHM 5 km s" 1 ). 

On every night, the calibration spectra were taken im- 
mediately before and after the scientific exposures. This en- 
sures that any environmental changes occurring during the 
observations are accurately accounted for. Since December 
2002, the UVES was equipped with an automatic reset- 
ting of the Cross Disperser encoder positions at the start 
of each exposureQ. The reason is to use the daytime ThAr 
calibration frames for saving the night time. This is fully 
justified for standard observations, but not for the measure- 
ments at the sub-pixel level which require the best possible 
wavelength calibration. To avoid the spectrograph resetting 
at the start of every exposure, the calibration frames were 
taken in the attach mode. The wavelength solutions were 
determined from each of the ThAr exposures and applied to 
the corresponding quasar exposures to make the pixel-to- 
wavelength conversion. For each frame of about 400 ThAr 
lines, more than 55% of the lines in the region were used to 
calibrate the lamp exposures. A polynome of the 5th order 
was adopted. Residuals of the wavelength calibrations were 
typically of 25 ms" 1 and symmetrically distributed around 
the final wavelength solution at all wavelengths. However, 
this represents the formal precision of the calibration curve 
but not the real calibration accuracy. 



1 Program ID 083.A.0733(A) 

2 http:/ /www. eso.org/observing/dfo/quality/UVES/pipeline/ 
pipe_reduc.html 
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Table 1. Journal of the observations. In the first column, calibration exposures are not numbered (blank). T m i n , T max are minimum 
and maximum air temperature during the exposure, and Pstart, fend are air pressure at the beginning and end of the exposure. 
ENC - encoder readings. 



No. of 


Date 


Starting 


Setting, 


Exposure, 




T 

max; 


-fstart j 


-fend? 


ENC 


cxpubuic 




lime, U J. 


A (nm) 


sec 


K 


K 


mmHg 


mmHg 






2009-09-20 


01:38:16 


437 


25 


12. 


.5 


12 


.5 


743 


.53 


743. 


53 


1079359 




2009-09-20 


01:38:12 


760 


1 


12. 


.4 


12 


.5 


743 


.53 


743. 


.53 


557957 


1 


2009-09-20 


01:39:46 


760 


5400 


12. 


.4 


12 


,5 


743 


.57 


743. 


.75 


557957 


1 


2009-09-20 


01:39:50 


437 


5400 


12. 


.4 


12 


.5 


743. 


.57 


743. 


.75 


1079359 




2009-09-20 


03:11:00 


437 


25 


12 


.4 


12 


,4 


743 


.73 


743. 


.72 


1079359 


2 


2009-09-20 


03:12:35 


760 


5400 


12. 


.4 


12 


,4 


743 


.77 


743. 


63 


557954 


2 


2009-09-20 


03:12:39 


437 


5400 


12. 


.4 


12. 


.4 


743 


.77 


743. 


.63 


1079359 




2009-09-20 


04:43:46 


437 


25 


12. 


.4 


12 


.4 


743 


.60 


743. 


.62 


1079359 




2009-09-20 


04:43:43 


760 


1 


1 

Iz. 


A 

.4 


1 

Iz 


A 

.4 


7 A Q 
1 4o 


.DU 


7/1 Q 
( 4o. 


Dz 


557954 




2009-09-20 


06:46:15 


437 


25 


12 


1 


12 


2 


743 


03 


743 


03 


1079357 


3 


2009-09-20 


06:47:36 


437 


5400 


11 


Q 

. tJ 


12 


1 


743 


00 

-VJVJ 


742 


70 


1079357 


3 


2009-09-20 


06:47:32 


760 


5400 


11 


q 


12 


1 


743 


00 


742 


70 


557947 




onnn no on 
zUUy-U9-zU 


Uo.l0.4o 


A Q7 
4o / 


zo 


12 


1 


12 


2 


749 


68 


742 


68 


1U (9oO ( 




2009-09-21 


01:44:27 


580 


25 


11. 


.7 


11 


.7 


744 


.13 


744. 


,18 


4690970 




2009-09-21 


01:44:30 


390 


25 


11. 


.7 


11 


.7 


744. 


.13 


744. 


.18 


1093943 


4 


2009-09-21 


01:46:01 


580 


5400 


11. 


.7 


11 


.7 


744 


.17 


743. 


.88 


4690970 


4 


2009-09-21 


01:46:05 


390 


5400 


11. 


.7 


11 


.7 


744 


.17 


743. 


.88 


1093943 




2009-09-21 


03:17:09 


580 


25 


11. 


.7 


11 


.7 


743. 


.88 


743. 


.88 


4690970 




2009-09-21 


03:17:12 


390 


25 


11. 


.7 


11 


.7 


743 


.88 


743. 


.88 


1093943 


5 


2009-09-21 


03:18:44 


580 


5400 


11. 


.0 


11 


7 
. 1 


7 A Q 
1 4o. 


.00 


7/1 Q 
( 4o. 


DO 


4690970 


5 


2009-09-21 


03:18:48 


390 


4507 


11 


6 


11 


7 


743 


88 


743 


68 


1093943 




2009-09-21 


04:35:12 


580 


25 


11 


6 


11 


6 


743 


70 


743 


68 


4690970 




2009-09-21 


04:35:14 


390 


25 


11 


6 


11 


.6 


743 


.70 


743 


68 


1093943 




2009-09-22 


00:12:33 


420 


25 


11. 


.7 


11 


.7 


742. 


.88 


742. 


.88 


1084636 


6 


2009-09-22 


00:14:11 


420 


5400 


11. 


.1 


11 


.7 


742 


.90 


743. 


.15 


1084636 


6 


2009-09-22 


00:14:07 


700 


5400 


11. 


.1 


11 


.7 


742 


.90 


743. 


.15 


544738 




2009-09-22 


01:45:20 


420 


25 


11. 


.1 


11 


.1 


743. 


.15 


743. 


.15 


1084636 


7 


2009-09-22 


01:46:58 


700 


5400 


10. 


.8 


11 


.1 


743 


.17 


742. 


.78 


544727 


7 


2009-09-22 


01:47:02 


420 


5400 


10. 


.8 


11 


.1 


743 


.17 


742. 


.78 


1084634 




2009-09-22 


03:18:12 


420 


25 


10. 


.8 


10. 


.8 


742 


.78 


742. 


.80 


1084634 


8 


2009-09-22 


03:19:55 


420 


4279 


10 


.8 


10 


.8 


742 


.78 


742 


.48 


1084632 


8 


2009-09-22 


03:19:51 


700 


4200 


10. 


.8 


10 


.8 


742 


.78 


742. 


.50 


544720 




2009-09-22 


04:32:19 


420 


25 


10. 


.8 


10 


.8 


742 


.48 


742. 


.50 


1084632 




2009-09-22 


04:32:10 


700 


1 


1U 


.6 


1U 


Q 


7 A O 


A Q 
.40 


7/1 
/ 4z 


A Q 
.40 


544720 




2009-09-23 


00-21-28 


580 


15 


10 


.4 


10. 


.4 


742 


.30 


742. 


.30 


4690939 




2009-09-23 


00:21:31 


390 


25 


10. 


.4 


10 


.4 


742. 


.30 


742. 


30 


1093936 


9 


2009-09-23 


00:22:51 


580 


5400 


10. 


,3 


10 


,4 


742. 


.33 


742. 


.38 


4690939 


9 


2009-09-23 


00:22:54 


390 


5400 


10. 


.3 


10 


.4 


742 


.33 


742. 


.38 


1093936 




2009-09-23 


02:03:37 


390 


25 


10 


,3 


10. 


,3 


742 


.38 


742. 


.40 


1093936 




2009-09-23 


02:03:34 


580 


15 


10 


.3 


10 


.3 


742. 


.38 


742. 


.40 


4690939 


10 


2009-09-23 


02:05:17 


580 


3600 


10 


.2 


10 


.3 


742 


.38 


742. 


.15 


4690933 


10 


2009-09-23 


02:05:18 


390 


3600 


10. 


.2 


10 


.3 


742 


.38 


742. 


.15 


1093935 




2009-09-23 


03:06:27 


580 


15 


10. 


.2 


10. 


.2 


742 


.17 


742. 


.20 


4690933 




2009-09-23 


03:06:29 


390 


25 


10. 


.2 


10 


.2 


742 


.17 


742 


.20 


1093935 


11 


2009-09-23 


03:07:50 


580 


3600 


10. 


,2 


10. 


.2 


742 


.13 


741. 


.87 


4690933 


11 


2009-09-23 


03:07:52 


390 


3864 


1U 


.2 


1 A 

10 


.2 


7 A O 

742 


.13 


7/11 

(41. 


.83 


1093935 




2009-09-23 


04:13:13 


580 


25 


10. 


.2 


10 


,2 


741. 


.83 


741. 


.87 


4690933 




2009-09-23 


04:13:23 


390 


25 


10. 


.2 


10 


.2 


741. 


.83 


741. 


.87 


1093935 


12 


2009-09-24 


02:14:01 


580 


3740 


11. 


.1 


11 


.3 


741 


.87 


741. 


.88 


4690957 


12 


2009-09-24 


02:14:05 


390 


3740 


11 


.1 


11 


,3 


741 


.90 


741. 


.88 


1093940 




2009-09-24 


03:17:29 


580 


15 


11. 


.1 


11 


.1 


741 


.90 


741. 


.88 


4690957 




2009-09-24 


03:17:33 


390 


15 


11. 


.1 


11 


.1 


741 


.90 


741. 


.88 


1093940 


13 


2009-09-24 


03:19:21 


580 


3740 


11. 


.1 


11 


.1 


741. 


.88 


741. 


.88 


4690954 


13 


2009-09-24 


03:19:25 


390 


3740 


11. 


.1 


11 


.1 


741 


.88 


741. 


.88 


1093939 




2009-09-24 


04:22:49 


580 


15 


11. 


,1 


11. 


.1 


741 


.90 


741. 


.90 


4690954 




2009-09-24 


04:22:52 


390 


25 


11. 


.1 


11 


.1 


741 


.90 


741. 


.90 


1093939 


14 


2009-09-24 


05:17:51 


580 


3740 


11. 


.0 


11 


.1 


741 


.70 


741. 


.50 


4690952 


14 


2009-09-24 


05:17:55 


390 


3740 


11. 


.0 


11 


.1 


741 


.70 


741. 


.50 


1093939 




2009-09-24 


06:21:19 


580 


15 


11. 


.1 


11 


.1 


741 


.50 


741. 


.48 


4690952 




2009-09-24 


06:21:22 


390 


15 


11. 


.0 


11 


.1 


741 


.50 


741. 


.48 


1093939 
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Table 2. Wavelength ranges covered by individual settings 



Setting AA, A 



390 


3290 


-4519 


420 


3586 


-4820 


437 


3759 


-4985 


5801 


4788 


- 5763 


580u 


5836 


- 6809 


7001 


5530 


- 6936 


700u 


7050 


- 8904 


7601 


5694 


- 7532 


760u 


7660 


- 9465 



The precision of the wavelength calibration in the course 
of science reduction is determined by the temperature and 
air pressure gradients between the times of the science ex- 
posure and of the wavelength calibration applied to the 
data since the environmental changes may cause drifts in 
the refractive index of air inside the spectrograph between 
the ThAr and quasar exposures. These thermal-pressure 
drifts move the different cross dispersers in different ways 
thus introducing the relative shifts between the different 
spectral ranges for different exposures. For instance, ac- 
cording to Kaufer et al. (2004) the pressure changes in the 
UVES enclosure induce a shift of w 66.7 m s^ 1 /mmHg, 
whereas temperature fluctuations produce shifts of w 420 
m s _1 /degr C in the red arm and w 60 m s _1 /degr C in 
the blue arrr0. The magnitudes of the possible shifts in our 
dataset can be estimated from the minimum and maximum 
temperatures inside the spectrograph enclosure during the 
exposures and from the pressure values at the beginning 
and the end of each exposure which are given in Table [T] 
In the last column of this table the encoder readings are 
reported which indicate that there were no grating resets 
within each block of observations. We emphasize that this 
effect has never been taken into account in all the studies 
performed so far on UVES data. There are no measurable 
temperature changes for the short exposures of the cali- 
bration lamps but during the science exposures the tem- 
perature drifted generally by 0.1-0.2 K. Pressure changes 
range from 0.15 to 0.6 mmHg (note that the start and the 
end pressure values are inverted in the UVES fits headers). 
Thus, measured temperature and pressure changes assure 
a radial velocity stability within 50 m s _1 in the blue 
arm and within w 100 m s _1 in the red arm. 

Individual spectra were corrected for the motion of the 
observatory around the barycenter of the Earth-Sun sys- 
tem and reduced to vacuum. The component of the obser- 
vatory's barycentric velocity in the direction to the object 
was calculated using the date and time of the integration 
midpoint. The air wavelengths were transformed to vacuum 
by means of the dispersion formula by Edlen (1966). We 
note that the barycentric and vacuum corrections require a 
rebinning which introduces a certain degree of correlations 
between the fluxes in the adjacent pixels. 

3. Analysis 

The laboratory line wavelengths used throughout the paper 
are listed in Table |3l For Mgi, Mgn, and Fen lines (except 

3 http:/ /www. eso.org/observing/dfo/quality/UVES/pipeline/ 
pipe_reduc.html 



Table 3. Atomic data for Mg I, Mg II, Al II, Al in, Si II, Ca I, 
Can, Fei, and Fell resonance transitions 



Line 


\ * A 

' v vac i 




f ^vac 


Mgi 


2852.96282(10) 1.83 


1 


Mgn 


2803.53112 


s 


0.306 


1 


Mgn 


2796.35403 


1(1 


0.615 


1 


Alii 


1670.7886 




1.74 


O 

z 


Aim 


1862.7910 




0.278 


2 
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References: [1] - Aldenius (2009); [2] - Morton (2003); 
[3] - Nave & Sansonetti (2010). 

*A V ac absolute errors affect the second to last digit in [2]. 
1 Oscillator strengths / are taken from [2]. 



for A1608 A), we adopted the wavelengths from Aldenius 
(2009). The wavelengths of the Mgi line from Salumbides 
et al. (2006), the Mgn lines from Batteiger et al. (2009), 
and the Fen lines from Nave & Sansonetti (2010) all are 
shorter by < 0.2 mA than the corresponding values from 
Aldenius (2009). In the velocity space, the shifts between 
these wavelength scales are in the range from 3.5 m s _1 
to 26 m s _1 , but since all wavelengths are shifted in the 
same direction, the relative shifts between lines which are 
compared in our analysis do not exceed 15 m s _1 which is 
well below the measurement errors. The photospheric solar 
abundances are taken from Lodders et al. (2009). 

The ThAr wavelength calibration of QSO spectra 
taken with slit echelle spectrographs such as, e.g., the 
Keck/HIRES or the VLT/UVES, is not stable and can 
produce both velocity offsets between different exposures 
and intra-order velocity distortions within each exposure 
(Levshakov et al. 2006; Centurion et al. 2009; Griest et al. 
2010; Whitmore et al. 2010). The reason of these calibra- 
tion distortions is still not well understood. They can be 
caused by differences in the slit illumination between the 
source and the calibration lamp, by non-linearity of echelle 
orders, non- uniformity in the distribution of the reference 
ThAr lines, and by other still obscure factors. In general, 
the detected velocity shifts are of order of a few hundreds 
of m s _1 , i.e. well below the pixel size. Thus, if the ef- 
fects to be studied have a characteristic scale of about the 
pixel size, QSO spectra can be prepared by means of com- 
mon pipeline procedures. In our case the pixel size is ~ 1.3 
km s _1 and the expected line shifts are about hundreds of 
m s^ 1 , i.e., we work on a sub-pixel level. In order to guar- 
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Fig. 1. An example of the cross-correlation of two profiles 
from different exposures. Upper panel: Normalized intensi- 
ties of the Mg II A2796 A line from the z a b s = 1-5864 system 
(Sect. I4.2[) extracted from exposure #2, setting 760 (solid 
line) and from exposure #1, setting 760 (dotted line). See 
Table Q] Lower panel: The dependence of S - the sum of 
squares of the intensity differences over all pixels covering 
the line profiles - as function of the velocity shift, Sv. The 
profiles become fully aligned when the dotted line profile is 
shifted by w —350 m s _1 (s 1/4 pixel size). 



antee the required accuracy, the profiles of all absorption 
lines used in the present study were prepared individually 
in accord with the following procedure. 

At first, we calculate the redshift z a b s of a given 
absorption-line system. This redshift is fixed for all other 
lines from this system. All line position measurements are 
carried out in the velocity scale. Namely, for each absorp- 
tion line with the laboratory wavelength Ao, a spectral re- 
gion centered at A bs = Ao (1 + z a t>s) and wide enough 
to incorporate clear 'continuum windows' from both sides 
(blueward and redward) of A b s is converted into the ve- 
locity scale as V — c(A — A bs)/A bs, where c is the speed 
of light and A is the current wavelength. The continuum 
windows are used to determine the local continuum level 
and to normalize the extracted line profiles. 

As already mentioned above, the profiles from individ- 
ual exposures may be shifted relative to each other. To eval- 
uate these shifts we used a procedure that cross-correlates 
two normalized profiles 1\ and I2 '■ one of the chosen profiles 
(Ii) is fixed and the other (I2) is moved relative to I\ by a 
small step, Sv (usually Sv w 50 m s _1 , i.e., of about 0.05 
pixel size). At each step, a value of S - the sum of squares 
of the intensity differences over n pixels covering the line - 



is calculated, S(Sv) = Y^=i[^,j ~ h,j(Sv)] 2 . S depends 
parabolically on Sv, being minimal when the profiles are 
aligned. This shift at the minimum of S is taken as a local 
velocity offset, AV, between the two profiles. An example 
of such cross-correlation is shown in Fig. [1] 

Cross-correlating the profiles of different lines from dif- 
ferent exposures we found that (i) the velocity shifts be- 
tween individual exposures of the same line were, in gen- 
eral, within |AV| < 500 m s _1 , (ii) the shift between two 
arbitrary chosen exposures was not a constant but varied 
from line to line (i.e., wavelength dependent). This supports 
the above statement that calibration accuracy is affected by 
many factors. Then, a natural solution to obtain the result- 
ing line profile would be to average over all factors, i.e., to 
sum up all available exposures. 

However, a particular feature of the present observa- 
tions is that because of the seeing conditions the signal 
changes significantly between exposures of different set- 
tings. Namely, for the setting 390+580, 8 exposures were 
taken, all with a comparable signal level. But for the setting 
437+760, only 3 exposures were obtained, each of which 
with a signal much higher than in the setting 390+580. 
The high-signal exposures have larger weights in the sum 
and, hence, even a single miscalibrated high-signal profile 
can noticeably shift the combined profile. 

To avoid such shifts we firstly co- added only profiles 
with a comparable signal and then checked whether the 
profiles from the high signal exposures - if present - were 
shifted with respect to this combined profile or not. An ex- 
ample is shown in Fig. [2j left panel. The shifted exposures 
were aligned with the combined profile and then all avail- 
able profiles were co-added together to provide the final 
line. 

Calculations showed that the alignment of the high sig- 
nal exposures resulted in a shift of the final line center of 
about 40-50 m s _1 as compared to a line obtained by a 
simple co-adding (i.e., without alignment) of all profiles. 
This value is of the same order as the statistical errors 
with which the centers of strong and narrow absorption 
lines can be measured (see, e.g., Tables 0HS]) and, there- 
fore, the alignment is important while dealing with such 
lines. On the other hand, the centers of weak lines are mea- 
sured with errors of 100-200 m s _1 which means that the 
alignment of the high signal exposures does not lead to sta- 
tistically significant differences. Nevertheless, in the present 
study the profiles of weak lines were also prepared with the 
alignment since in many cases this operation allowed us to 
obtain more regular (smooth) profiles. However, the pro- 
files of weak lines can be severely distorted by noise which 
would result in artificially exaggerated AV values (Fig. [2 
right panel). That is why the exposure shifts were calcu- 
lated for a strong absorption line located in the vicinity 
(within 15 — 20 A) of a weak line and then applied to the 
individual profiles of this weak line. 

The wavelength rang e above 7200 A is covered by the 
3 exposures of the setting 760 (divided into low and upper 
frames) and by the 3 exposures of the setting 700 (also low 
and upper frames), all having significantly different signals. 
For absorption lines from this range, the individual profiles 
were aligned with a profile from the exposure with the high- 
est signal. In turn, the local calibration (i.e., at the position 
of a particular line) of this exposure was checked using dif- 
ferent transitions of the same ion (see Sect. 14.21 and 14.4( 1. 
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Fig. 2. Cross-correlations of profiles from exposures with high signal levels. All lines are from the z a b s = 0.45207 system. 
Left upper panel: Solid line - the normalized Mg II A2796 A profile composed from 12 exposures with similar signal levels 
(100— 150 counts per pixel at the continuum level); dotted line - the normalized Mgn A2796 A profile from the exposure 
#2, setting 437 (Tabled]) with a high level of signal (~ 500 counts per pixel). Left lower panel: S as function of Sv (same 
as in Fig.[lJ. The dotted line profile is to be shifted by « —250 m s _1 to become fully aligned with the composite profile. 
Right upper panel: Solid line - the normalized profile of the weak Cai A4227 A line composed from 8 exposures with 
similar signal (~ 300 counts per pixel); dotted line - the normalized Cai A4227 A profile from the exposure #1, setting 
760 (Table [1]) with a high level of signal (~ 1000 counts per pixel). Right lower panel: S as function of Sv (same as in 
Fig- HI)- The dotted line profile is distorted by noise which leads to its apparent shift relative to the composite profile of 
about —1000 m s _1 . The real shift of this portion of the exposure #1 is only —150 m s _1 as calculated from profiles of 
a strong Fen A2382 A line at z a b s = 1.5864 located in the vicinity of Cai A4227 A. 



The combined profiles were normalized by means of a 
smoothing spline and then the rms noise was calculated us- 
ing the continuum windows. These final profiles were fitted 
to Gaussians defined by the standard parameters: the com- 
ponent center V, the broadening Doppler parameter 6, and 
the column density N. The FWHM values to convolve the 
synthetic profiles with the spectrograph point-spread func- 
tion were estimated from the ThAr lines in the vicinity of 
particular absorption lines. These values range between 4.3 
and 5.4 km s _1 . However, in some cases a simultaneous fit- 
ting of several lines of the same ion shows that the FWHM 
should be slightly (~ 5%) smaller than the value estimated 
from the ThAr lines. This can be explained by additional 
broadening of the ThAr lines caused by the pressure effects 
since these lines arise in dense plasma. 

The errors of the fitting parameters were calculated via 
the inverse second derivatives matrix. Since we are espe- 
cially interested in the accurate error of the line center V, 
this error was additionally evaluated using the Ax 2 -method 
(Press et al. 1992). As the actual error of the line center the 
larger estimate was accepted. 



It is well known that the noise produced by CCDs is 
highly correlated (e.g., Levshakov et al. 2002). These cor- 
relations do not affect the results of the ^-minimization, 
but they reduce the apparent dispersion of the noise and, 
hence, lead to the underestimation of the calculated errors 
of the fitting parameters. In order to understand how to 
correct the estimated errors, the following Monte Carlo ex- 
periment was performed. First, we synthesized a line profile 
with a given set of parameters. The spectrum of HE0001- 
2340 reveals wide continuum windows (e.g., in the range 
4474-4692 A or 5130-5314 A) which in fact represent long 
sequences (several thousands of points) of correlated noise 
values. In these windows, we choose at random a point and 
starting from it extracted n subsequent values which were 
added to the synthesized profile (n is the number of pixels 
covering the line). Then the profile parameters were evalu- 
ated via the standard ^-minimization. The procedure was 
repeated 1000 times, and the dispersions of the resulting 
distributions of the fitting parameters were taken as their 
errors. They were w 1.5 times larger than the errors esti- 
mated via the inverse matrix (or the errors estimated from 
the Ax 2 -procedure with the confidence level of A% 2 = 1). 
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Fig. 3. The observed profiles of the Fe II and Mg II 
lines from the z a b s = 0.45207 system towards the quasar 
HE0001-2340. The system consists of two subcomponents 
separated by « 70 km s _1 . The zero radial velocity is fixed 
at z = 0.45207. Note the difference between relative line 
intensities in the subcomponents. 

Thus, the final errors listed in Tables 2HH1 below represent 
the calculated values multiplied by 1.5. 

The difference between the radial velocities V\ , Vi of 
two absorption lines can be used to estimate a putative 
variation of the fine-structure constant a in space and time 
with respect to its laboratory (terrestrial) value. For two 
transitions with different sensitivities to a changes the value 
of Aa/a = {a z — ai a b)/ai a b is calculated as Aa/a = (Vi — 
V2)/[2c(Q2 — Qi)}, where c is the speed of light, and Qi, Q2 
are the corresponding sensitivity coefficients (Dzuba et al. 
2002; Berengut et al. 2010). 



4. Results 

4.1. Absorber at z a h s 



0.45207 



This unique systems exhibiting a set of Fei lines was dis- 
covered and firstly described by D'Odorico (2007). The sys- 
tem consists of two components separated by « 70 km s _1 
(Fig. [3]). The red component shows only a strong Mgn 
doublet and weak Fen AA2600, 2382 A lines, but the blue 
component contains, along with Fei, narrow and well dis- 
tinguished lines of Mgn AA2796,2803 A, Mgi A2852 A, 
Fen AA2600, 2382, 2344, 2374 A, Can AA3969,3934 A, and 
Cai A4227 A with simple profiles and represents in fact the 
best-choice system to measure the Mg isotope abundance 
ratio. In the previous spectrum of HE0001-2340 taken with 
a lower spectral resolution (pixel size ~ 2.3 km s _1 , FWHM 
= 6.8 km s" 1 ), the Mn triplet AA2576, 2594, 2606 A can be 
safely identified and probably also a line of Sil A2515 A. 
Unfortunately, in our spectrum which was obtained with a 
higher resolution and because of this has a lower S /N these 
lines are indistinguishable from noise fluctuations. 

Calculations were performed with the line profiles pre- 
pared as described in Sect. [3] Some special cases are out- 
lined below. 

Fei A2967 A (Fig. @|. The line is weak and not seen in 
individual exposures. However, this Fe 1 transition has the 




Fei 2484 

J_i 



S/N = 24 





Fei 2523 



S/N = 22 



Fei 2967 



S/N = 5E 



Fei 3021 



S/N = 58 



_FeI 3720 



S/N 



40 



-40 -30 -20 -10 



Relative radial velocity (km s _1 ) 

Fig. 4. Fei lines from the z a b s = 0.45207 system towards 
the quasar HE0001-2340 (solid-line histograms). Synthetic 
profiles are plotted by the smooth curves. The signal-to- 
noise ratio, S/N, per pixel in the co-added spectrum is in- 
dicated in each panel. The zero radial velocity is fixed at 
z = 0.45207. 



highest sensitivity to changes in a (Dzuba & Flambaum 
2008), so it was important to obtain its line profile as accu- 
rate as possible. Fortunately, 12 A redward to the position 
of Fe 1 A2967 (at A « 4309 A) lies a strong Al 11 A1670 A line 
from the z a b s = 1.5864 system (Sect. S3) which can be used 
to calculate the shifts of the individual exposures and to test 
the local wavelength scale calibration. The corresponding 
spectral ranges are present in all 14 exposures - the 8 expo- 
sures of the setting 390 with similar signal, the 3 high-signal 
exposures of the setting 437, and the 3 intermediate-signal 
exposures of the setting 420. Employing our standard pro- 
cedure of the line profile preparation, we would sum up 
the 8 equal-signal profiles from the setting 390, align the 
profiles from the high-signal exposures with this combined 
profile, and then co-add all available exposures. I.e., the 
reference for the final profile would be set by the average 
profile from the setting 390. But the comparison of the Al II 
A1670 A line produced by the 8 co-added exposures of this 
setting with other low-ionization lines Si 11 A1526 Aand Fe II 
AA2382, 2600, 2586, 2344, 2374, 1608 A from the same z abs = 
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Fig. 5. Same as Fig. @] but for Fell lines. The right panel shows the decomposition of the A2586 A line which is blended 
with two Ly-a forest absorption lines (see text, for detail). 

1.5864 system showed that the Alii profile was shifted rela- 
tive to these lines by i=s 0.5 km s _1 , whereas the profile com- 
posed from the combined 6 exposures of the settings 437 
and 420 was coherent with them. Therefore, in the present 
case, the reference profile was obtained by co-adding the ex- 
posures from settings 437 and 420, and the alignment was 
performed for the individual Al n profiles from the setting 
390. The same procedure was used to prepare the resulting 
profile of the Fe I A2967 A line: the exposures of the settings 
437 and 420 remained unchanged, whereas those from the 
setting 390 were shifted in accord with the shifts calculated 
for the Al II A1670 A profiles. Additionally, exposures ex- 
hibiting noise spikes within the Fei A2967 A profile were 
excluded from the final co-adding (3 exposures from the 
setting 390 and 3 from the setting 420). 

Fei A3021 A (Fig.QJ. In spite of its weakness, the line is 
well distinguishable in all available 14 exposures (settings 
390, 437, 420). However, from unknown reasons (may be 
because of the absence of appropriate ThAr lines in the 
vicinity of A = 4387 A where Fei A3021 is located) the 
final profile appeared to be offset relative to the other Fe I 
lines. However, the sensitivities to a variations of the Fe I 
A3021 A and A2967 A transitions are equal, and these lines 



should be aligned at any rate. The cross-correlation of the 
Fe I A3021 A and A2967 A lines showed that the former was 
shifted by -0.35 km s" 1 . The calibration of the Fei A2967 
A line was verified by the Alii A1670 A line from the z a b s 
= 1.5867 system, and therefore we aligned the Fei A3021 
A line with Fei A2967 A, i.e., shifted it by 0.35 km s _1 . 

Fei A3720 A (Fig. g]). The line is weak and present 
only in the 8 exposures of the frame 5801. Because of the 
similar sensitivity to a changes, its position should coin- 
cide with the Fe I A2484 A and Fe I A2523 A lines, but the 
cross-correlation showed that the Fe I A3720 A line was off- 
set by —0.3 km s" 1 . On the other hand, the profiles of 
Fei AA2484, 2523, 2967, and 3021 (corrected) were coherent 
which means that Fe I A3720 was affected by some miscal- 
ibration. Thus, the final line profile of Fei A3720 A was 
shifted by 0.3 km s _1 . 

Fen A2586 A (Fig ]E\ right panel). This quite strong 
line is blended with Ly-a forest absorptions, but can be 
accurately deconvolved. The noise was calculated as the 
difference between the observed and the fitted spectrum 
and added to the deconvolved Fe n A2586 A profile shown 
in Fig. [SJ left panel. 
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Fig. 6. Same as Fig. H but for Mgl, Mgll, Cal, and Call lines. 

Table 4. Fitting parameters of spectral lines identified in the z a b s = 0.45207 absorber. Indicated arc la statistical errors. 



Ion 


V, km s- 1 


b, km s 1 


N, cm~ 2 


C 


Mgn 


-2.21 ±0.04 


2.79 ± 0.25 


(2.0 ± 0.1)E12 


0.77 


Mgl 


-1.94 ±0.17 


2.43 ± 0.39 


(1.9 ± 0.1)E11 


0.77* 


Fen 


-1.77 ±0.06 


2.59 ± 0.05 


(1.06±0.03)E13 


0.86 


Fcl 


-1.72 ±0.14 


0.84 ±0.10 


(2.9 ± 0.2)E12 


0.71 


Can 


-1.82 ±0.16 


1.52 ± 0.19 


(6.7±0.4)E11 


0.56 


Cai 


-1.84 ±0.27 


1.6 ± 1.0 


(7.2 ± 1.8)E10 


0.56* 



'Single line, covering factor C adopted from the correspondin 
ionized doublet. 



Cai A4227 A (Fig. EJ). This weak line is located at A = 
6138 A. Close to its position, there are a weak Fen A2374 
A (at A = 6142 A), and a strong Fe n A2382 A (at A = 6163 
A) lines, both from the z a bs= 1-5864 system. The lines are 
present in the 14 exposures: 8 from the frame 580u, 3 from 
7001, all 11 with a similar signal, and 3 from the frame 
7601 with high but different in all 3 exposures signals. The 
velocity shifts of the profiles from the frame 7601 relative to 
the co- added 11 profiles from the frames 580u and 7201 were 
calculated for the high-contrast Fen A2382 A line. Then 
the same shifts were applied to the Fe II A2374 A exposures 
from the setting 7601. The final profiles of Fe n AA2382, 2374 
A were perfectly aligned with the Fell AA2600, 2586, 1608 
A lines and with other low-ionization lines from the z a b s 
= 1.5864 system. The same procedure was applied to the 
Cai line: profiles from the frame 7601 were shifted in accord 
with the Fen A2382 A exposures and then co-added with 
unchanged profiles from the frames 580u and 7001. 

Mgl A2852 A (Fig. [5]). The line is weak and located at 
4142 A. Close to its position at 4150 A, there is a high- 



contrast line of Oi A1032 A from the z a b s = 2.1871 system 
(Sect. l4~4)) . The Mg I profiles from the high-signal exposures 
(settings 437) were corrected following the shifts of this O I 
line. 

Mgn AA2796, 2803 A (Fig.©. Both lines are strong and 
present in all 14 exposures. However, the profiles of Mgn 
A2803 A in 4 exposures (2 of the setting 390, 1 of the setting 
437, and 1 of the setting 420) are severely distorted and not 
included in the final spectrum. This explains different S/N 
ratios for the Mgn A2796 A and Mgn A2803 A lines shown 
in Fig. El 

As a result, we have for the analysis the follow- 
ing set of the absorption lines: Mgl A2852, Mgn 
AA2796,2803, Fei AA2484, 2523, 2967, 3021, 3720, Fen 
AA2600, 2382, 2586, 2344, 2374, Cai A4227, and Can 
AA3934, 3969 A. The line profiles of each ion were analyzed 
independently. Trial calculations have shown that a single- 
component Gaussian profile is sufficient to fit properly all 
lines. However, the value of x 2 < 1 in the fitting of the 
Mgn doublet (highest S/N ratio) could be achieved only 
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Fig. 7. The velocity shift AV as function of the ratio £ = ( 25 Mg+ 26 Mg)/ 24 Mg for the Mgi A2852 A (upper panel) and 
the Mgn A2796 A (lower panel) lines. The curves are calculated for different ratios of the heavy isotopes, r = 26 Mg/ 25 Mg 
(depicted over each curve). Ay = corresponds to the solar value of the Mg isotope ratios, £ Q = 0.2660 and r wl. 
The horizontal dashed lines mark the mean values of AV MgI and AV MeU measured at z a b s — 0.45207, whereas the dotted 
lines restrict the ±ler uncertainty intervals (see text, for detail). 



with a covering factor C < 1, i.e., assuming an incomplete 
coverage of the background light source: alternative 
models with two components and C — 1 yield very small 
6-parametcr (« 0.5 km/s) which is inconsistent with the 
line widths of other species from this system. Spectra of 
other ions have lower S/N ratios and their fits with C < 1 
and with C = 1 are statistically indistinguishable, although 
X 2 in trials with C < 1 was always lower (e.g., 0.7 vs. 0.9). 
Thus, the final calculations were performed with C < 1 for 
all ions. 



The covering factor C is easily determined from 
a slightly modified model of the line profile, I\ = 
Cexp(— r\) + (1 — C), under the assumption that C is the 
same for the whole absorbing cloud (i.e., it does not de- 
pend on A or, equivalently, on the radial velocity V). Here, 
/ is the observed normalized intensity within the absorp- 
tion line profile and t\ is the optical depth. The covering 
factor can differ from ion to ion since, in general, different 
ions trace different gas. If at least two lines of the same ion 
are available, then the covering factor is estimated safely. 
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Fig. 8. Empirical probability distribution functions (p.d.f.) 
for Mg isotopic fractional abundances, /, estimated by a 
Monte Carlo procedure (see Sect. 14. 1[ for detail) 



If only one line is present, then the covering factor cannot 
be determined (degenerate with column density) and has 
to be fixed to some value (see note to Table S]). In fact, the 
incomplete covering is quite common among the so-called 
associated (i.e. located close to the background quasar) ab- 
sorbers and not such a rare finding even among intergalactic 
absorbers (see examples in Paper I) . The cause of the light 
leakage is not clear, but, probably, it is related to the grav- 
itational lensing of the quasar beam by some intervening 
galaxy/cluster. In any case we can conclude that the ab- 
sorber showing the incomplete covering of the background 
source is to be very small, probably even of a sub-parsec 
linear size. Note that the incomplete covering in the z a b s = 
0.45207 system was recently mentioned also by Jones et al. 
(2010). 

The observed and synthetic profiles from the z a b s = 
0.45207 system are shown in Figs. an d the correspond- 
ing fitting parameters are listed in Table |U 

The radial velocity of the system center calculated as 
the weighted mean value of the Fell, Fel, Call, and Cai 
line centers is V c — —1.77 ± 0.03 km s _1 . Then for the 
velocity offsets of the Mg lines we obtain the following 
values: AV MgII 

^^Me I = I 



Ktgn - Vc. = -0.44 ± 0.05 km s" 1 



V c = -0.17 ± 0.17 km s" 1 , i.e., nega- 
tive velocity shifts are detected for both the Mg n lines and 
the Mgi line although for the latter it is only at the la- 
level because of a large uncertainty of the Mgi line posi- 
tion (weak single line). However, before ascribing this shift 
to the enhanced content of the heavy Mg isotopes in the ab- 
sorbing gas, the arguments concerning possible calibration 
errors and kinematic shifts (the so-called 'Doppler noise') 
should be considered. 

The centers of both Fe n and Fe I lines are determined 
from 5 transitions, and therefore they provide a robust 
mean radial velocity. The calibration of Ca II and Ca I lines 
is verified by the Fe n AA2382, 2374 A lines from the z a b s = 
1.5864 system and their centers coincide with those of the 
Fen and Fel lines within the la uncertainty interval. 



Both Mg II lines fall in the spectral region with the high- 
est S/N ratio. They are narrow and strong - this explains 
a small error (0.04 km s^ 1 ) of the line center determina- 
tion. However, the systematic error due to miscalibration 
can easily reach the tenfold value as was shown above for, 
e.g., the Fei A3021 A and Fel A3720 A lines. Can the posi- 
tions of the Mg II lines be also affected by some additional 
errors ? The answer is very probable negative. The main 
factor here is that the position of Mgn A2803 coincides 
with a strong and unsaturated line Thl A4070 A from the 
ThAr spectrum which is used as a reference for the wave- 
length calibration. As a result of this coincidence, the Mgn 
A2803 profiles from different exposures are almost aligned 
(shifts < 0.05 km s _1 ). Individual profiles of the Mgn 
A2796 reveal shifts up to ±250 m s" 1 (see, e.g., Fig. H]), 
but the final Mgll A2796 profile prepared by the standard 
procedure (i.e., similar-signal exposures co-added to form a 
preliminary profile, high-signal exposures aligned with this 
profile and then all exposures co-added again) is fully co- 
herent with Mgn A2803. In general, all lines in the present 
study which were found to be shifted due to calibration 
errors (e.g., Fei A3021, this Section, Fen A1608 at z a b s = 
2.1871, Sect. 14. 4[) fall in the gaps between ThAr reference 
lines, whereas lines with good calibration (e.g., Sin A1526 
and Fell A2344 from the z abs = 2.1871 system, Sect. l4~4l) 
coincide with some of the ThAr lines. The conclusion that 
the vicinity to some of the ThAr reference lines ensures a 
stable wavelength calibration is further supported by the 
comparison of the line positions in the archive HE0001- 
2340 spectrum and in the present one: the 'coinciding' lines 
are coherent in both spectra whereas centers of lines from 
the gaps can differ by ~ 0.5 km s _1 . Thus, the negative 
shift of the Mg n lines in the system under consideration 
cannot be caused by calibration errors. The same statement 
is valid for the Mgi A2852 line as well, since its position is 
also close to a ThAr reference. 

Neither can this shift be attributed to the gas flows 
inside the absorber. The estimated Doppler ^-parameters 
(Table 2]) show that the line broadening is not thermal and, 
hence, some turbulent component is present. But the line 
centers of Fe II and Fe I (as well as Ca n and Ca i) coincide 
within the la uncertainty interval - contrary to the centers 
of the Mg II and Mg I lines. However, at values of the ioniza- 
tion parameter log U < — 4, which is the range most prob- 
able for the system under study (D'Odorico 2007; Jones 
et al. 2010), the ionization curves of the Mgn and Fell 
ions are almost parallel whereas the ionization curves of 
Fei and Fen diverge. Thus, in case of any density/velocity 
gradients inside the absorbers one would expect rather a 
shift between the centers of the Fe I and Fe n lines but not 
between Fen and Mgn. 

It remains to assume that the measured negative shift 
of the Mgll and Mgi lines relative to other lines in the 
Zabs = 0.45207 system is indeed due to enhanced content 
of heavy Mg isotopes in the absorbing gas. Figure [7] shows 
the velocity shift AV as a function of the ratio £ = ( 25 Mg+ 
26 Mg)/ 24 Mgfor the Mgi A2852 A (upper panel) and the 
Mgn A2796 A (lower panel) lines. The curves are calculated 
for different ratios of the heavy isotopes, r — 26 Mg/ 25 Mg, 
which are depicted over each curve. The laboratory wave- 
lengths for separate isotopes are taken from Salumbides et 
al. (2006) and from Batteiger et al. (2009) for, respectively, 
Mgi and Mgn. AV — corresponds to the solar value 
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Fig. 9. The overview of the H I Ly-a, Fe II A2600 A, Mg II A2796 A, C IV A1548 A, and Si IV A1393 A profiles from the 
£abs= 1-5865 absorption complex towards the quasar HE0001-2340. The complex consists of four sub-systems (^4, B, C, D) 
spread over 1100 km s^ 1 . The zero radial velocity is fixed at za = 1.5864. 



of the Mg isotope ratios, £© = 0.2660 and r© w 1. The 
horizontal dashed lines mark the measured mean values of 
AVmbi = -170 m s" 1 and AV Ms u = -440 m s^ 1 , and the 
dotted lines restrict the ±lcr uncertainty interval for Mgn 
(for Mgi, it runs from AV MgI — to —340 m s" 1 ). From the 
lower panel of Fig. [7] we see immediately that the measured 
shift of the Mg n lines can be produced by the isotope mix- 
ture with £ > 1.5 and r > 1, i.e., with an overabundance of 
26 Mg relative to other two isotopes. 

Using the values of AV MgI and AV MgU and their er- 
rors, we can estimate the most probable fractional abun- 
dances /24,/25j and /26 of Mg isotopes assuming that the 
measured velocity shifts are random values normally dis- 
tributed with the means — 170 m s _1 and —440 m s _1 and 
the dispersions of 170 m s _1 and 50 m s _1 , respectively. 
Fig.[5]shows the empirical probability distribution functions 
for each isotope obtained by Monte Carlo modeling when 
AVJtfgi and AV M&11 were chosen at random from the cor- 
responding normal distributions and /24,/25> and /26 were 
calculated under the conditions /24 + /2s + /26 = 1> each 
fk > 0, £ = (/ 2 5 + /26V/24 > 1-5, and r = / 26 // 2 5 > 1- It 
is seen that the abundances of the first two isotopes 24 Mg 
and 25 Mg are almost rectangularly distributed, whereas the 
fractional abundance of 26 Mg is distributed normally. The 
calculated statistics are /24 = 0.19±0.11, /2s = 0.22±0.13, 
and / 26 = 0.59 ± 0.06. 

In order to understand what causes the enhanced abun- 
dances of heavy Mg isotope in the z a b s = 0.45207 system the 
origin of this absorber should be clarified. Consider now the 
column densities of the ions presented in Table |U Although 



Table 5. Fitting parameters of spectral lines identified in 
the z a bs = 1-5864 absorber. For each ion, the upper row 
lists parameters of the main component whereas parameters 
of auxiliary components are listed below. Indicated are la 
statistical errors. 



Ion 


V, km s" 1 


b, km s 


N, cm" 2 


Fen 


3.56 ±0.05 


1.92 ±0.04 


(3.48 ± 0.03)E12 




-6.75 ± 1.26 


4.5 ± 1.0 


(1.28 ± 0.12)E11 




9.59 ±0.30 


5.2 ± 1.6 


(1.34±0.12)E11 


Sin 


3.52 ±0.05 


2.09 ±0.10 


(1.03±0.01)E13 




-5.2 ±5.5 


19.0 ± 1.5 


(8.69 ± 1.161E11 




5.7 ±0.3 


8.8 ±0.1 


(6.41 ± 0.15)E12 


Alii 


3.55 ±0.15 


3.05 ±0.30 


(1.15 ± 0.05)E12 




-12.0 ±2.1 


14.1 ± 1.2 


(1.3 ± 0.1)E11 




7.4 ± 1.1 


7.4 ± 1.0 


(1.78±0.18)E11 


Aim 


3.65 ±0.36 


3.24 ±0.30 


(1.19 ± 0.04)E12 




10.41 ± 1.51 


11.57 ±2.02 


(5.28 ± O.lOjEll 




-11.08 ±5.32 


22.48 ± 5.20 


(1.02 ± 0.06)E12 


Mgn 


3.71 ±0.04 


2.79 ± 0.02 


(9.10 ± 0.10)E12 




-5.81 ±0.53 


21.13 ± 0.66 


(1.43±0.04)E12 




5.23 ±0.10 


9.65 ±0.11 


(2.25 ± 0.03)E12 


Mgi 


3.58 ±0.45 


3.74 ±0.82 


(6.1±0.6)E10 



accurate ionization corrections are not known, we can un- 
ambiguously conclude that the relative abundance of iron 
to magnesium, (Fe/Mg) > 5, is by almost an order of 
magnitude higher than its solar value, (Fe/Mg)© = 0.83 
since the ionization corrections for Fe 11 are always higher 
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Fig. 10. Metal absorption lines from the z a bs= 1-5864 system towards the quasar HE0001-2340 (solid-line histograms). 
Synthetic profiles are plotted by the smooth curves. The signal-to- noise ratio, S/N, per pixel in the co-added spectrum 
is indicated in each panel. The zero radial velocity is fixed at z = 1.5864. 



than that for Mg II irrespectively of the shape of the ion- 
izing radiation and the value of the ionization parameter 
U . The ionization curves of magnesium and calcium at 
l°g U < — 4 closely trace each other, and from the compar- 
ison of the column densities of Mg I and Mg II with those 
of Ca I and Ca II we obtain that the relative content of cal- 
cium to magnesium, (Ca/Mg) ~ 0.4, considerably exceeds 
its solar value, (Ca/Mg)© = 0.06. As already mentioned 
above, in the system under study there is probably the Si I 
A2515 A line, but its identification is uncertain since it falls 
into the Ly-a forest. If the absorption at the expected po- 
sition of the Si I line is indeed due to this transition, then 



iV(Sii) = (7.0 ± 0.7) x 10 11 cm" 2 (D'Odorico 2007). The 
ionization correction for Si I at log U < — 4 is very close 
to that for Mgl, and comparing iV(Sii) with JV(Mgi) we 
obtain again a significant relative overabundance of sili- 
con to magnesium: (Si/Mg) = 3.7 vs. (Si/Mg) Q = 0.98. 
Such an abundance pattern - a high iron content, overabun- 
dances of Si and Ca relative to Mg - is usually observed in 
remnants of type la supernovae (e.g., Stehle et al. 2005; 
Mazzali et al. 2008; Tanaka et al. 2010). The Mnii triplet, 
although indistinguishable in our data, is clearly seen in the 
low-resolution spectrum of HE0001-2340 where N(Mn n) 
= (3.6 ± 0.4) x 10 11 cm -2 was determined by D'Odorico 
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Table 6. Fitting parameters of spectral lines identified in 
the z a bs = f -6515 absorber. Indicated are la statistical er- 
rors. 



Ion 


v, km s 1 


b, km s 


N, cm -2 


Sill 


-4.00 ±0.30 


4.0 ± 1.4 


(1.75 ±0.14)E12 


Fen 


-4.10 ±0.30 


2.2 ± 1.0 


(2.35±0.50)E12 


Mgn 


-4.13 ±0.03 


2.53 ±0.03 


(4.00 ± 0.04)E12 


Mgi 


-4.33 ±0.39 


2.7 ± 1.5 


(6.46 ± 1.42)E10 



(2007) . It is not clear how to translate this column den- 
sity into the Mn abundance, but taking into account the 
trace solar content of Mn, (Mn/Fe)© = 0.01, a very high 
relative abundance and even an overabundance of Mn rel- 
ative to Fe seems to be quite probable. And this supposes 
a high metallicity progenitor, since, in general, the abun- 
dance of Mn in the supernova la remnants is thought to 
be correlated with the metallicity (expressed as the initial 
content of 22 Ne) of the exploding white dwarf (McWilliam 
et al. 2003; Badenes et al. 2008). Model calculations of nu- 
cleosynthetic yields from SN la do not predict enhanced 
content of heavy Mg isotopes for the whole SN la remnant, 
but very high concentration of 26 Mg is predicted for the 
outer layer of the exploding white dwarf where it is syn- 
thesized via alpha-capture on 22 Ne (Seitenzahl et al. 2010). 
The SN remnants are known to be heterogeneous objects 
displaying complex variations of chemical composition with 
position and velocity (e.g., Hughes et al. 2000). In fact, the 
z a bs = 0.45207 system also consists of two separated compo- 
nents (Fig. [3]) with obviously different physical parameters. 
Taking into account a tiny linear size (manifested through 
the incomplete coverage of the light source) of the analyzed 
clump at V — km s _1 we can suggest that it is a fragment 
of a SN la remnant. 

As already mentioned above, different sensitivity of dif- 
ferent Fei, Fen, and Can transitions to changes in the fine- 
structure constant a makes it possible to estimate Aa/a 
from the velocity shifts of the corresponding line centers. 
Unfortunately, the measured errors of the line centers in 
the present system allows us to put constraints on Aa/a 
only at the level of 10 -5 , with the most accurate estimate 
provided by Fei lines: Aa/a — (0.6 ± 1.0) x 10~ 5 . 

We note that Chand et al. (2004) and Murphy et al. 

(2008) estimated Aa/a in this system using the Mgn 
AA2796,2803 A and Mgi A2852 A transitions as the ref- 
erence. A negative Aa/a value reported by Murphy et al. 
is obviously a consequence of the unaccounted Mg isotope 
shift. 

4.2. Absorber at z ahs = 1.5864 

This system represents a large absorbing complex extend- 
ing over - 1000 km s" 1 (Fig. [9]). The redshift z = 1.5864 
corresponds to the sub-system A with the strongest low- 
ionization lines (Fen, Sin, Mgn, Alii). The sub-system 
D at V = —1200 km s _1 was described in Paper I as 
a high-metallicity absorber enriched by AGB-stars. The 
sub-systems A and B probably originate in low-ionization 
clumps (lines of Fn, Mgn, Sin, Cn, Alii) embedded in 
gas with a higher ionization which is seen in lines of Si IV 
and Civ. Note a remarkable kinematic similarity between 
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Fig. 11. Metal absorption lines from the z a bs= 1-6515 
system towards the quasar HE000 1-2340 (solid- line his- 
tograms). Synthetic profiles are plotted by the smooth 
curves. The signal-to-noise ratio, S/N, per pixel in the co- 
added spectrum is indicated in each panel. The zero radial 
velocity is fixed at z = 1.6515. 



Mgn and Fen absorptions in these sub-systems and Mgn 
and Fe II absorptions in the previously described system at 
z abs = 0.45207 (see Fig.©. 

Here we describe the sub-system A in the context of 
measuring the Mg isotope abundance ratio. Unfortunately, 
in the sub-system B the Mgn AA2796,2803 lines are con- 
taminated with telluric absorptions which prevent to mea- 
sure accurately their centers. The observed line profiles 
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shown by histograms in Fig. 1101 were prepared as described 
in Sect. [3l Some particular details are given below. 

Fell A2586 A. Three profiles from the frame 580u are 
severely corrupted and were excluded from the final co- 
added spectrum. This explains a lower S/N at the position 
of this line as compared to the neig hboring Fell A2600 A 
line. 

Fe II A2600 A. The cross-correlation of the Fe II A2600 A 
and Fe II A2382 A profiles revealed a shift of 0.15 km s _1 be- 
tween them. Since trial fittings of Fe II lines have shown the 
consistency of Fen A2382 A with Fen AA2344, 2586, 2374, 
and 1608 A, the final Fen A2600 A line was aligned with 
Fen A2382 A (i.e., shifted by 0.15 km s" 1 ). 

Fen A1608 A. The line is very weak and barely visible 
in the individual exposures. In the spectrum of HE0001- 
2340, it lies at 4160 A - close to a strong Oi A1302 A line 
from the z abs = 2.1871 system at 4150 A (Sect. S3}- The 
position of this O I line is verified by comparison with the 
Sin A1260 A and Sin A1526 A lines from the z abs = 2.1871 
system. Thus, the Fell A1608 A profiles from 3 exposures 
of the high-signal setting 437 were shifted in accord with 
the shifts of the O I A1302 A line. 

Alii A1670 A. See Sect. Q where Fei A2967 A line is 
described. 

Mg II AA2796, 2803 A. Both lines are very strong and 
present in 3 exposures of the frame 7601 and 3 exposures 
of the frame 700u, all having different S/N ratios. The ex- 
posure with the highest S/N (#8 from 700u, Tabled]) was 
taken as a reference and all other exposures were aligned 
with it before co-adding. 

Mg i A2852 A. The line is very weak and its profiles in 
the individual exposures (3 from 7601 and 3 from 700u) 
are severely distorted by noise. All available exposures 
were simply co-added without any preliminary alignments. 
This means that the final profile is mostly affected by the 
strongest exposure #8 from 700u. 

C II A1334 A. The line is almost saturated, has a low 
S/N and is partly blended with Ly-a forest absorptions. 
The line center cannot be estimated with accuracy better 
than 1 km s^ 1 and we did not include the C II A1334 A line 
in the final analysis. 

The lines of each ion (Fen, Sin, Alii, Aim, Mgn, Mgi) 
were fitted independently using Gaussian components. In 
each case the number of components was determined as a 
minimum required to fit the profile with % 2 < 1. The fitting 
parameters are listed in Table [5j the synthetic profiles are 
shown in Fig. [TU]by the solid lines. The weighted value for 
the system center calculated from the Fen, Sin, and Alii 
lines is V c = 3.54 ± 0.03 km s _1 . This gives the offset for 
the Mgn line of AF Mg „ = V UgU -V c = 0.17 ±0.05 km s" 1 
and AV MgI = 0.04 ± 0.45 km s" 1 for Mgi. 

As mentioned above, the profiles of the Mg n lines from 
individual exposures were aligned with the exposure #8 
from the frame 700u having the strongest signal. The offset 
revealed is the sum of two effects: a shift of the reference 
exposure, and a shift of the Mgn lines due to enhanced con- 
tent of heavy isotopes. The value of the exposure shift can 
be estimated from the following consideration. The Mgn 
lines are located in the range AA = 7230 - 7252 A . The 
frame 700u starts at 7050 A (Table [2J , and up to the po- 
sitions of the Mgn lines there are no absorption lines in 
the spectrum which can be used to verify the wavelength 
calibration. However, the frame 7601 (another frame whose 
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Fig. 12. Absorption lines from the z a b s = 2.1871 sys- 
tem overlapped to demonstrate different profile shapes. The 
Fen A2600 A line is shifted by 0.35 km s^ 1 , and the Mgn 
A2796 A line by 1.08 km s _1 with respect to their original 
positions in the spectrum (see Sect. 14.41 for detail). The 
zero radial velocity is fixed at z = 2.1871. 

exposures contribute to the Mg II profiles) starts at shorter 
wavelengths (covering AA = 5694 — 7532 A), and up to 
6809 A it overlaps with the wavelength interval covered by 
the frame 580u. In the range covered by both the 7601 and 
580u frames there are several absorption lines from three 
different systems (Can AA3934,3969 A, Cai A4227 A at 
z abs = 0.45207; Fen AA2344, 2374, 2382, 2586, 2600 A from 
the present system; and Fen A2382 A at z a b s = 1-6515) 
which can be used to cross-check the calibration of the 
exposures from the frame 7601. This cross-checking shows 
that the strongest exposure #3 has a systematic positive 
shift from 0.3 km s" 1 (positions of Can AA3934, 3969 A, 
Cai A4227 A at z ahs = 0.45207) to 0.2 km s" 1 (positions 
of Fen AA2586,2600 A from the present system) relative 
to the center of each absorption-line system. In turn, the 
Mgn profiles from the exposure #8 (700u) and #3 (7601) 
coincide, i.e., at the position of the Mgn lines these two 
exposures are coherent. Thus, we can assume that the ref- 
erence exposure #8 (700u) has a positive shift of < 0.2 
kms" 1 . This gives us a conservative lower limit on the 
offset of the Mg II lines due to enhanced content of heavy 
isotopes of AV MsII > -0.08 km s" 1 . Assuming 26 Mg/ 25 Mg 
> 1, which is supported both by measurements in presolar 
spinel grains and in some Milky Way giants (Gyngard et 
al. 2010; Yong et al. 2006) and by theoretical predictions 
(Karakas et al. 2006), we obtain ( 25 Mg+ 26 Mg)/ 24 Mg < 0.7 
(see Fig. [71 lower panel). We note that in the previous (lower 
resolution) spectrum of HE0001-2340 the Mgn lines from 
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Fig. 13. Metal absorption lines from the z a b s = 2.1871 system towards the quasar HE0001-2340 (soliddine histograms). 
Synthetic profiles are plotted by the smooth curves. The signal-to-noise ratio, S/N, per pixel in the co-added spectrum is 
indicated in each panel. Also shown are the relative velocity offsets AV: for the Sill A1190 A line and the Mgll doublet 
they are with respect to the Sin AA1260, 1526 A lines, and for the Fen lines - relative to the Fen A2344 A line. The zero 
radial velocity is fixed at z = 2.1871. 



the z a bs = 1.5864 system are shifted with respect to the 
Mgll lines from the present spectrum by w —0.5 km s" 1 . 

Again, the presence of spectral lines with different sen- 
sitivities to changes in a makes it possible to estimate 
the value of Aa/a, but now with a considerably higher 
accuracy since the errors of the line centers are small. 
Namely, by comparing the Sin A1526 A line as a tran- 
sition unaffected by a-variations with a sensitivity coef- 



ficient Q w 0.001 (Dzuba ct al. 2002) with the Fen 
AA2600, 2586, 2382, 2374, 2344 A transitions, all having sim- 
ilar sensitivity coefficient of Q w 0.04 (Porsev et al. 2007), 
we obtain Aa/a = (—1.5 ± 2.6) x 10~ 6 . This is compara- 
ble with the most stringent constraint on \Aa/a\ < 2 ppm 
found from the Fell lines at z = 1.15 towards the bright 
quasar HE 0515-4414 (Quast et al. 2004; Levshakov et al. 
2005, 2006; Molaro et al. 2008b). The line Fen A1608 A 
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with Q = —0.016 is very weak and, besides, its red wing 
is slightly corrupted by the noise spike (Fig. [TO]) . If only 
the blue (incorrupted) portion of the Fe n A1608 A profile 
is taken, then Aa/a estimated from the Fen transitions is 
Aa/a = (-0.5 ± 11.0) x 10" 6 . 



Table 7. Fitting parameters of spectral lines identified in 
the z a bs = 2.1871 absorber. For each ion, the upper row 
lists parameters of the main component whereas parameters 
of auxiliary components are listed below. Indicated are la 
statistical errors. 



4.3. Absorber at z a b s = 1.6515 

This system was described in Paper I where the following 
relative abundances were determined: [C/H] — 0.4 — 0.45, 
[Mg/C] = 0.1 - 0.15, [Fe/C] < - 0.3, and [Si/C] = -0.8. 
Such an abundance pattern relates the absorbing gas with 
the outflows from AGB-stars. A high content of magnesium 
along with a relative deficit of iron is one of the characteris- 
tics of hot-bottom burning (Herwig 2005) which can result 
in an enhanced content of the neutron-rich Mg isotopes. 
Indeed, analyzing the previous (lower resolution) spectrum 
of HE000 1-2340 we noticed that Mgn lines were offset by 
w —0.6 km s _1 relative to the lines of other low-ionization 
species from this system. Here we present accurate line cen- 
ters based on the higher resolution spectrum. 

The available low-ionization transitions are Cn A1334 
A, Sin AA1260, 1526 A, Fen A2382 A, Mgn AA2796,2803 
A, and Mgi A2852 A. However, the Cn A1334 A and Sin 
A1260 A lines are blended with strong Ly-a absorptions and 
cannot be deconvolved with a sufficiently high accuracy. 
The combined Mgn AA2796, 2803 A profiles consist of 6 
co-added exposures (3 from 7601, and 3 from 700u) aligned 
with the exposure #8 (700u). The weak Mgi A2852 A line 
is present only in 3 exposures of the frame 700u which were 
simply co-added, i.e., the profile of this line is mostly af- 
fected by the exposure #8 with the strongest signal. The 
profiles of Sill A1526 A and Fen A2382 A were prepared as 
described in Sect. |3 

The line centers were determined for the Sin A1526 A, 
Fe ii A2382 A, Mg n AA2796, 2803 A, and Mg i A2852 A lines. 
Lines of each ion were fitted independently. The blue wing 
of Mg II A2796 A is blended with a weak telluric absorption, 
but the central part of this line is unaffected. The fitting 
parameters are listed in Table [6l the observed and synthetic 
profiles are shown in Fig. [TT] The system center calculated 
as the weighted mean between the centers of the Si n A1526 
A and Fen A2382 A lines is V c = -4.05 ± 0.20 km s" 1 
which gives the offsets AV Mgll = -0.08 ± 0.20 km s" 1 , and 
AV Mg i = -0-28 ± 0.44 km s" 1 . 

The Mgil lines are observed in the range AA = 7414 — 
7434 A. At both sides of this region there are absorption 
lines with verified calibration. Namely, at 7380 A, the Mgi 
A2852 A line from the z a b s = 1.5864 system is present 
and its position coincides with other lines from this sys- 
tem (Table IU). At 7472 A, we have the Fen A2344 A line 
from the z a b s = 2.1871 system, and this line is also un- 
shifted (see Sect. I4.4j) . Thus, the calibration of the Mgn 
lines is correct. As for the Mgi A2852 A line, it lies close 
to the Fen A2374 A line from the z abs = 2.1871 system 
which is shifted relative to Fen A2344 A by 0.2 km s _1 
(Sect. [14]). Therefore, the Mgi A2852 A line is very proba- 
bly shifted as well, i.e., V MeI should be — 4.13±0.39 km s _1 , 
and AV Ug i = —0.08 ± 0.44 km s _1 . Because of weakness 
of the Sill A1526 A and Fen A2382 A lines their positions 
are measured with large errors making the shift of Mg lines 
statistically indistinguishable from zero. Under assumption 



Ion 


V, km s- 1 


b, km s 1 


N, cm -2 


Fen* 


4.94 ±0.10 


2.56 ±0.05 


(6.12±0.04)E12 




2.54 ±0.15 


6.47 ±0.17 


(2.09±0.05)E12 




13.73 ±0.35 


2.88 ±0.29 


(4.06±0.14)E11 


Sin 


4.92 ± 0.09 


2.58 ±0.06 


(9.92±0.35)E12 




-6.6 ± 1.0 


9.13 ±0.47 


(5.94±0.17)E11 




3.15 ±0.15 


6.57 ±0.05 


8.77±0.13)E12 




14.83 ± 0.10 


4.55 ±0.07 


(1.2±0.1)E12 


Oi 


4.97 ±0.05 


4.49 ± 0.03 


(3.10±0.13)E14 




-6.7 ±3.1 


1.0 ± 1.0 


(1.73±0.05)E12 




18.9 ±0.5 


4.8 ±2.2 


(2.85±0.65)E12 


Mgn 


3.84 ±0.10 


3.75 ±0.06 


(7.76±0.14)E12 




-7.14 ±0.33 


10.88 ±0.41 


(7.10±0.26)E11 




1.15 ± 0.15 


6.97 ±0.05 


(4.56 ±0.05)E12 




13.67 ±0.15 


4.35 ±0.11 


(8.22±0.10)E11 


'Lines 


are aligned with respect to Fen A2344 A (see Fig. [T3J). 



that 26 Mg/ 25 Mg > 1, a tentative upper limit on ( 25 Mg 
26 Mg)/ 24 Mgis < 2.6. 



4.4. Absorber at z ahs = 2.1871 

This absorber belongs to an extended sub-DLA system de- 
scribed in detail in Richter et al. (2005). It exhibits multiple 
lines of O I, Si II, Mg II, and Fe n which can be used to check 
the calibration uncertainties. 

The Oi AA1302, 1039 A, Sin AA1526, 1304, 1260, 1190 
A, Fen AA2600, 2586, 2382, 2374, 2344, 1608 A, and Mgn 
AA2803, 2796 A lines are present in many exposures of dif- 
ferent settings. Their profiles were prepared as described in 
Sect. 03 

The Fen AA2586,2600 A lines are detected in 3 expo- 
sures of the frame 760u and in 3 exposures of the frame 
700u (all exposures with different S/N). Before co-adding, 
the exposures were aligned with the exposure #3 (760u) 
showing the highest S/N. We note that Fell profiles from 
this exposure completely coincide with Fen profiles from 
the exposure #8 (700u) which is the highest S/N exposure 
in the setting 700u. 

The Fe n A2344 A line is present in 3 exposures of 7601 
and in 3 exposures of 700u. Again, all of them have different 
S/N. The alignment was performed relative to the exposure 
#8 (700u). 

The Fe n AA2382, 2374 A lines are present only in 3 ex- 
posures of 700u where their profiles are completely con- 
sistent, and the final profiles were obtained by simple co- 
adding. 

The Mgn AA2796, 2803 A lines are present in 3 expo- 
sures of 760u. Before co-adding, individual profiles were 
aligned with the profile from the exposure #3. 

Line profiles of different ions in the z a b s = 2.1871 sys- 
tem differ from each other which is clearly seen in Fig. [12] 
where several profiles are overplotted. This means that the 
absorbing gas has density and velocity gradients. However, 
all ions trace some central condensation, and we can esti- 
mate the shifts between the lines of different ions using their 
central components. Lines of each ion were fitted indepen- 
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dently to Gaussians the number of which was determined as 
a minimum required to fit all transitions of a given ion with 
X 2 < 1. The synthetic profiles are shown by solid lines in 
Fig. [T3] The fitting parameters are listed in Table Some 
details on the fitting procedure are explained below. 

The ion O I is represented by the lines 1302 A and 1039 
A. The O I A1302 A line is almost saturated and the 1039 
A line falls in a very noisy part of the spectrum. To con- 
strain the fitting parameters we included in the analysis the 
O I A988 A line from the previously obtained spectrum of 
HE0001-2340 with FWHM = 6.8 km s" 1 (its profile is also 
shown in Fig. I13[) . 

Trial fittings of Sin lines showed that the Sin A1190 
A line was noticeably shifted with respect to the Sin 
AA1260, 1526, 1304 A lines and that the synthetic Si n A1304 
A line came out systematically weaker than the observed 
one (see Fig. [13]). Both the 1190 A and 1304 A lines were 
excluded from the fitting procedure, and the fitting param- 
eters were estimated on base of the 1260 A and 1526 A 
lines. After the fitting procedure, the synthetic profile of 
the 1190 A line was calculated and the shift of the observed 
1190 A line was evaluated by cross-correlation of the syn- 
thetic and observed profiles. We found that Sin A1190 A 
line was shifted by —0.6 km s _1 relative to the other Sin 
transitions. This line is located just at the order edge what 
probably explains the revealed discrepancy. As for the Si n 
A1304 A line, the reason for its intensity to be underesti- 
mated whereas the profiles of the Sin A1260, 1190 (shifted), 
and 1526 A lines are fitted perfectly, is unclear. May be it 
is due to a blend with some unidentified absorption. 

Trial fittings of the available Fen lines 
(AA2600, 2586, 2374, 2382, 2344, 1608 A) revealed that 
all of them are shifted with respect to each other. We 
are especially interested in calibration of the Fen A2344 
A line since it lies in the vicinity of the Mgn doublet 
from the z a b s = 1.6515 system (Sect. I4.3|) and can be 
used to verify its position. We took Fen A2344 A as a 
reference line, calculated offsets for the other Fen lines, 
and aligned all profiles. Then the aligned Fen lines were 
fitted together. The fitting parameters and the shifts of 
the individual lines relative to Fen A2344 A are shown in 
Table [7J The center of the 2344 A line coincides within the 
uncertainty interval with the centers of the O I and Si n 
lines. Thus, we can conclude that the exposure #8 from 
the frame 700u (reference exposure for the alignment of 
the individual Fen A2344 A profiles before co-adding) at 
7472 A (position of Fen A2344 A) has no offsets relative 
to the zero-point. Note that there is a good ThAr reference 
line at 7472 A. However, negative offsets begin to appear 
at larger wavelengths being —0.2 km s _1 at 7564 A (Fen 
A2374 A), -0.25 km s" 1 at 7594 A (Fen A2382 A), -0.45 
km s- 1 at 8244 A (Fen A2586 A), and -0.35 km s" 1 at 
8286 A (Fen A2600 A). 

A negative offset is clearly displayed also by the Mgn 
AA2796, 2803 A doublet located in the range 8912-8936 A: 
it is shifted by —1.08 ± 0.13 km s _1 relative to Sin lines 
(Table [7J). To what extent this offset is due to a putative 
high content of heavy Mg isotopes is not clear. According 
to Richter et al. (2005), in the present system [O] = —1.8 
and [N/O] < —1.5, i.e., enrichment occurred due to metal- 
poor SN II. Thus, from general considerations the enhanced 
ratio ( 25 Mg+ 26 Mg)/ 24 Mg should not be expected. It is 



likely that the shift of the Mg n lines is mostly caused by 
the distortions in the calibration process. 

This tendency - increasing negative offsets with increas- 
ing wavelength in the range A > 7500 A - is detected not 
only in the present spectrum, but in other UVES spectra as 
well. Such offsets are found both in the archive QSO spectra 
and in spectra newly obtained within the framework of the 
ESO Large Program 185.A-0745. At the hardware level the 
UVES is realized by the two arms - blue and red, and in the 
red one there are two CCD chips - lower (REDL) and upper 
(REDU). The blue arm covers the range 3000 - 5000 A, and 
the red arm 4200 -11000 A. The performance of the blue 
and lower red chips was checked by an original procedure 
involving solar radiation reflected by asteroids (Molaro et 
al. 2008a) and no large systematic shifts in the wavelength 
calibration between two chips were detected. Our data sup- 
port this result, although they also demonstrate that local 
line position shifts with amplitudes in most cases below 
|AU| = 0.4 km s _1 (< 1/3 pixel size) are possible. Why 
the wavelength calibration slides down in spectra detected 
with the REDU chip (7500 A < A < 9000 A) is unclear 
and this problem needs its thorough investigation. In any 
case it is quite obvious that the studies involving differential 
measurements of the line positions may use REDU spectra 
only if there are some independent methods to check their 
wavelength calibration. 

This is just the case of the system considered here 
where calibration of Fen A2344 A is verified by compar- 
ison with multiple lines of other ions (O I and Si n) . Since 
the sensitivity coefficients of O I is Q 1 « (Berengut 
& Flambaum 2010) and of Sin is Q SlII « 0.001 (Dzuba 
et al. 2002), positions of the Oi and Sill lines should 
not be affected by putative variations of the fine-structure 
constant a, whereas the sensitivity coefficient of the Fell 
A2344 A transition is 0.036 (Porsev et al. 2007). The dif- 
ference AV between the Fen A2344 A and Oi/Sin line 
centers is much smaller than the error of AV which is 0.11 
km s^ 1 if the Oi line center is taken as a zero-point, and 
0.13 km s _1 if Sin is a reference line. The former case 
gives us the value Aa/a = (1.4 ± 5.1) x 10~ 6 , the lat- 
ter - Aa/a = (-0.9 ± 6.0) x 10 -6 . In turn, the line Fen 
A1608 A with Q = -0.016 is shifted by 0.3 km s" 1 rel- 
ative to the Fen A2344 A, which would give us the value 
Aa/a = (9.6 ±4.5) x 10~ 6 . This example confirms the con- 
clusion of Griest et al. (2010) that calibration errors are 
the main source of uncertainties in Aa/a estimations from 
absorption lines in quasar spectra. 

We note that the present system was used by Chand et 
al. (2004) and by Murphy et al. (2008) to estimate Aa/a 
from the transitions Mgn A2803 A, Aln A1670 A, Sin 
A1526 A, and Fe n AA2374, 2586 A which were fitted simul- 
taneously using the same model. It is unclear how this could 
be done taking into account that in the previous HE0001- 
2340 spectrum all peculiarities described above (e.g., ex- 
treme shift of the Mgn lines and different line profiles of 
different ions) were also present. 

5. Conclusions 

In the present paper we attempt to measure the Mg isotope 
abundances in several absorption systems detected in the 
spectrum of the quasar HE0001-2340. The spectrum was 
obtained with the VLT/UVES with the slit width of 0.7 
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arcsec and was read pixel by pixel without binning, thus 
ensuring about two times smaller pixel size as was achieved 
in previous (archived) spectra of this source (1.3 km s _1 
vs. 2.3 km s" 1 ). The line profiles were prepared individu- 
ally according to a special procedure described in Sect. [3] 
Absorption-line systems selected for the study reveal mul- 
tiple lines of the same ions (e.g., Sin AA1526, 1304, 1260 A, 
Fen AA2382, 2374, 2344 A etc.) which can be used to verify 
the local wavelength calibration. The main results are the 
following. 

1. For the first time we measured the abundance of Mg 
isotopes in the high-redshift absorber at z a b s = 0.45207 
which is probably a remnant of the SN la explosion of 
high-mctallicity white dwarf (s) . Lines of the Mg II dou- 
blet are shifted relative to low-ionization transitions of 
other ions (Fei, Fen, Cai, Can) by AV UsU = — 0.44± 
0.05 km s" 1 , and the Mgi line by AV MsI = -0.17±0.17 
km s~ 4 . This translates into the isotope abundance ra- 
tio 24 Mg: 25 Mg: 26 Mg = (19 ± 11) : (22 ± 13) : (59 ± 6) 
with strong relative overabundance of heavy Mg iso- 
topes 25 Mg+ 26 Mg relative to 24 Mg as compared to the 
solar ratio 24 Mg: 25 Mg: 26 Mg = 79:10:11. 

2. In the absorbers at z a b s = 1.5864 and z a b s = 1.6515 
arising in the gas enriched likely by outflows from AGB 
type stars we obtain for the shift of the Mg n lines rela- 
tive to other ions the values of AV MgU > — 0.08 km s _1 
at z abs = 1.5864, and AV Msll = -0.08 ±0.20 km s" 1 at 
z&bs — 1.6515. This allows us to set only upper limits 
on the content of heavy Mg isotopes in the absorbing 
gas: ( 25 Mg+ 26 Mg)/ 24 Mg < 0.7 at z abs = 1.5864, and 
(25 Mg+ 26^/24^ < 2 6 at Zabs = x 6515 

3. In the present spectrum of HE0001-2340 the errors of 
the wavelength calibration produce velocity shifts of ab- 
sorption lines with amplitudes in most cases below 0.4 
km s _1 (1/3 pixel size). This statement is valid for the 
wavelength ranges covered by the blue and REDL chips 
of the UVES. However, absorption lines taken with the 
REDU chip are systematically shifted toward shorter 
wavelengths. The discrepancy increases with increasing 
wavelength and may reach —1 km s _1 (« 0.8 pixel size) 
in the wavelength interval 7500 A < A < 9000 A. 

4. In the absorption system at z a b s = 1.5864 with the ver- 
ified calibration of the metal absorption lines, we set a 
limit on the variation of the fine-structure constant at 
the level of Aa/a = (—1.5 ± 2.6) x 10~ 6 which is one 
of the most stringent estimates of this value obtained 
from optical spectra of QSOs. 

5. Comparison of statistical errors of the line position mea- 
surements with systematic errors due to miscalibration 
of the wavelength scale shows that the systematic error 
is the main source of uncertainties in the measurements 
of Aa/a from quasar absorption-line systems. 
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